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Abstract. We have random number of independent diffusion processes with absorption 
on boundaries in some region at initial time t = 0. The initial numbers and positions of 
processes in region is defined by Poisson random measure. It is required to estimate of 
number of the unabsorbed processes for the fixed time r > 0. The Poisson random measure 
depends on t and r — > oo. 



Consider the set of independent random diffusion processes ^k{t), k = 1,N, 
t > 0, (,kiQ) = Xk, Xk G Q C R^. We wish to investigate of distribution of the number 
of the processes £,k{t) which was into Q for all moments of time t < t. 

Let domain Q C i?'' be open connected region and it is limited by smooth surface dQ. 
All processes £,k{t) are diffusion processes with absorption on the boundary dQ. These 
processes are solutions of the following stochastic differential equations in Q 

d 

1=1 

x), a{t, x) : R+x R'^ ^ R'^. 
with an initial condition: ^(0) = Xk G D. 

Here the W^''\t) = {wf\t), l<i <d), 1 < fc < A are independent in totality d- 
dimensional Wiener processes. 

Thus, these processes have the identical diffusion matrices and shift vectors , but they 
have different initial states. 

Let Q is bounded and boundary dQ is Lyapunov surface C^^''*'-'. 

The initial number and positions of processes are defined by the random Poisson 
measure /i(-,T) in Q: 

P(MAt) = fc) - ^^^^^e-"(^-^' 

where to(-, r) is finitely additive positive mesure on Q for fixed r. 

This task was offered in [f] as the mathematical model of practice problem. The 
authores in article [2] investigated case when initial number and positions of diffusion 
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processes are defined by determinate limited measure N{B,t). Where the N{B,t) is 
equal to number of points Xk in a set B and N = N{Q, r) < oo for fixed r > 0. 
We consider the following case 

a{t,x) = a = {0, . . . ,0), bi{t,x) = bi = {bii,...,bid), l<i<d; 

d 

We define matrix a = B^B, B = (bij), 1 < i,j < d 

o = (cy), i < i,j < d and differential operator A : ^ ^vdSox'- 

l<i,j<d ' ' 

Let (7 be a matrix with the following property 

cTijZiZj > 

l<i,j<d 

Here /x, there is fixed positive number, and z= {zi, - ■■ ,Zci) there is an arbitrary real 
vector. 

This operator acts in the following space 

Ha = {u:u€ L2{Q) nAu€ L2{Q) n u{dQ) = 0} 

with inner product {u,v)a = (Am, w). Here (, ) is inner product in L2{Q). The opera- 
tor A is positive operator. 

It is known [3] that the following eigenvalues problem 

Au ^ -\u, u{dQ) = 
has infinity set of real eigenvalues Aj ^ oo and 

< Ai < A2 < • • • < As < • • • . 

The corresponding eigenfunctions 

/ll ) • • • ) /ini ) ■ ■ ■ ) /si ) • • • ) fsug 1' ' ' 

form complete system of functions both in Ha and 

^2(0) ■= {u : u & L2{Q) n u{dQ) = 0}. Here the number nk is equal to multiplicity 
of eigenvalue A^. 

We denote by r]{T) the number of remaining processes in the region D at time instant r. 

We also assume that cr-additivc measure i' is given on the E,y- algebra sets of 

Q, v{Q) < 00. All eigenfunctions fij-Q^ and all measures to(-, r) are (Si/, Sr) 
measurable. Here is system of Borel sets of B}. Let denotes the weak convergence 
of random values. 

Put 

ff(T) =exp(-^Ai) . 
Theorem 1. We suppose that m{-,T) satisfies the condition 

lim m(B, T)g(T) = v{B), B G T.^. 

T — i-OO 

Then rj{T) ^ rj i/ t 00 where rj has Poisson distribution function with the 

"1 

parameter a = J F(x)dv and F(x) = ^ fii{x)cii, cu = J fii{x)dx. 
Q j=l Q 
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Proof. We consider the following initial-boundary problem 

du 1 X ^ d'^u 

l<i,J<d ■' 

u{0, a;) = 1 if x € Q; 



u{t, x)=0 if xedQ, t>0 (2) 

It is known [4], that u(r, x) is equal to probability of remaining in the region Q at time 
instant r of a diffusion process from (1) which occurs at the point (0,a:) at the initial 
moment ( ^(0) = x, x & Q). We designate through = {x\, ■ ■ ■ , a;^) the initial position 
of fc-th process. We define the value of ^(t, 7fe). 

We define a particular solution of (2) in form 

u{t, x) = Ul{t)u2ix). 

The ordinary argumentaion leads to definition of joined constant A: 

2—^^ — _ _x 
Ui dt U2 

We obtain the following system of tasks due the latter one 

Au2 = -\u2; U2{dQ) = 0. (3) 



^ = -^t.i; u,{0) = l (4), 

It is clear that ui{t, A) = exp(— |A) is solution of (5) . The soluton of (3) was described 
above. We assume that system of functions {fij{x),i > 1,1 < j < rii} is orthonormalized 
with respect to space L2{Q). 

The general solution of problem (2) has the following form 

oo ^ n, 

U{t,x) = ^exp(--Aj) ^ Cjmfjmix) 
j=l m=l 

where coefficients Cjm are equal to coefficients of decomposition of initial value (unit) 
by system of functions fjm- Cjm — J fjm{x)dx. The Parseval - Steklov equality is true 

Q 

for these coefficients: 

oo "j 

J2T.^"m = \Q\- (5) 

m 

Put -^(a;) = J2 cufim{x)- The function F{x) is continuous and bounded function on 

_ m=l 

the Q. Since u{t, x) is probability,it is not diffucult to show that F{x) > for all x Q Q. 

Let M = sup F{x). We introduce the following sets 
xeQ 
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Bk,n = {xeQ: < F{x) < — ^ ^} 

n n 

Here < k < n — 1 and n > 1. 

Let us denote by Ck,n{T), 1 < k < n the number of unabsorbed processes at time 
instant r which occur in the region Bk,n at the initial time. These values are independent 
in totality by assumption. The distribution function of Cfe.n {t) is defined by the following 
formula 

oo 

PiCkA^) = = ^P(M(-Bfe,„,T) = d)x 
d=l 

I d 

X n"(^'^'^) n (i-w(^>7ij), «=o,i,.... 

l<ii,--- ,ii<d,im7^»j,TO7^j fe=l s=(+l,»s^(»ir" M)tim¥^ij 

Here a;,^. e Sfe,„. 
We set 

afe,n(r) = min u(r, a;), afc,„(r) = 1 - afe,ri(r); 

?'fc,n(T) = max M(T,a;), 6fe,„(T) = 1 - 6fe,„(T). 



Now 



d=l 

< P{Ck,n{r) = l)< 

oo df ID \ 

E jj'"' exp(-m(i3fe,„,r))g>i,„(r)a^-'(r) =:4,„a,r). (6) 



Further 



(m(Bfe,„,T)afc,„(r))' ^ (&fc,„(T)m(Bfc,„,T))'^ ' 
= exp(-m(Bfe,„,r))2_^ (d^^ " 

exp(-6fe,„(T)m(Bfe,„,r)); 



{d-l)\ 

u—i 

{■m{Bk,n,r)ak,n{T)y 



By analogy: 

{m{Bk,n,T)bk,n{T)y / \ /D ^^ /'7^ 

Ik,n{l'T) = — — ^, ex.p{-ak,n{T)m{Bk,n,T)). (7) 

We introduce the following generating functions 

l>0 
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Vfc,n(r, s) = J2 s^P{Ck,n{^) = 1), k = 0,n-l, < s < 1. 

By the construction, r}{T) can be represented in the form r}{T) = Ci,nH 1" Cn- !,«(''")• 

Thus 

n-l 

^{t,s) =J\_^Pk,n{T,s). (8) 
fe=0 

Combining (6)- (8), we conclude that 

exp{(safe,„(r) - 6fe,„(T))m(Sfe,„, r)} < 'Pk.n{r,s) < 
< exp{(s6fe,„(T) - ak,n{T))m{Bk,ri,T)} 

and 



exp < ^(sofe,„(T) - bk,n{T))m{Bk,n,r) > < iy?(r,s) < 



. /c=0 

n-l 



< exp{^(s6fc,„(T) - afe,„(r))m(Bfe,n,T)}. (9) 



Since function u(t, x) is continuous function in a; G Q, there exit a points a;*, a;* S -B^,^ 
such that the following equalities have place 



ak,n{T) = Mi(T,Ai)F(a;*) + ^■ui(r, Afc) ^ ckmfkm{x*), 



k>2 m=l 



bk,n{T) =Ml(T,Ai)F(a;*) + ^Ui(r,Afe) ^ Ckmfkm{x*), 



k>2 m=l 



here a;* := a;*(fc,n,T), a;* := x*{k,n,T). 

Now, we can rewrite the sums in exponentes from (9) in the following forms 



n-l 

r , 



J2 {sF{x,) - F{x*)) exp(--Ai)m(Bfc,„, r)+ 

fe=0 

nj 

Xl^^P(~i^i)'^(^'='"''^)Xl'^^P (^i^^J ^ X] Cjm(s/jm(a;*) -/jm(a;*)), (10) 

■ ~ j>2 m=l 

n-l 

^(sF(a;*) - F{x^)) exp(-^Ai)m(Sfe,„, r)+ 



;s=0 j>2 m=l 

n-l 



/c=0 



n— i "'j 

+ ^ exp(-^Ai)rn(Bfc,„, r) ^exp (-^(Aj - Ai)) ^ Cjm(s/j,„(a;*) - fjm{x»)), (11) 

/s=0 j>2 m=l 

We calculate limit of (10) if r — > oo. The first sum of (10) convergers to the following 
limit under the condition of theorem 
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n-1 



k=0 



fe=0 



This is difference of two integral sums which has the following limit under n ^ oo (see 



[5]) 



(s- 1) J F{x)v{dx). 
Q 



Put 



St(.x) = ^exp (-^(Aj - Ai)^ E Ckmfkm{x). 
j>2 m=l 

We consider sums of eigenfunctions in the form 



e(x,A)= E fU^) 



\i<X 



The following result is proved in monography [6,Thm. 17.5.3] 



sup e{x, A) < CA2 . 
xeQ 



Asymptotic characteristic of eigenvalues Xj under j — > oo is defined by the following 
inequalities [3, sec. 18] 



cij < Aj < C2 j , where ci , C2 = const. 

The latter one, (5) and Caushy-Bunyakovskii inequality lead to the following conver- 
gence under r ^ oo 



\sAx)\ < Eexp (-^(A, - Ai)) . E ^'m\ E /: 

j>2 \ m=l \ m=l 



<x/cEAjexp(--(A,-A)), 

j>2 



\ E '^im - 

^ TO=1 



V j>2 



j>2 m=l 



;2 ^0. 



Thus the second sum from(lO) convergences to zero. 

The similar considerations apply to (11). Proof is complete. 

Example. Now we apply the general approach to the particular case. 
We consider the case if Q is circle Q = {{x,y) : x^ + y"^ < ^ol- We assume that the 
diffusion processes occurs at the point {xk,yk) S Q at the initial time. 

The processes are described in Q by the following stochastic differential equations 
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d^{t) = J2hdwi{t) (12) 

i=l 

C(o) = ^0 = {xk,yk)- 

where bi = {a, 0), 62 = (0, cr) and 

W{t) = (wi{t), i ^ 1,2) is 2-dimensional Wiener process. 

We assume that the equation (12) defines a diffusion process with absorption on the 
boundary dQ = {{x, y, z) : + = r^}. 

In follows that the Jo {x) , Ji {x) are Bessel functions zero and first order. They are 
defined as the solutions of the following equations 

(fv 1 dy , 

— - H ^ + (1 ) = 0, 

dx^ xdx x^ ' 

y{xQ) = 0, (a;o = \/Ar); |y(0)| < 00; 

under n = and n = 1. 

The value of is equal to m- th root of the equation Jo(/u) = [7,8]. 
Let mes(-) denotes the Lebesgue measure. 
We set 



/(r) := exp 



We suppose that m(-,r) holds the condition 

m{-,T)f{T) ^ mes{-) if t — » 00. 
In this case the system of tasks (3), (4) has the following form 

Au2 = -IJ.U2, {x, y) e C; U2{x, y) = if x'^ +y^ = (13), 

^ = -^^u„ «i(0) = l. (14) 

According to general approach for construction of solution u{t, x, y) (see,for example, 
[7, sec. IV ] we rewrite the task of (13) in polar coordinates: U3(r, (p) := U2(r cos ip, r sin tp). 
The U3 is solution the following problem 

(9^M3 1 du3 1 d'^Us 
or r or r o<p 

^3(^0, </?) = 0. 



We obtain 



u{t,x,y) = u{t,r) = ^ CmJo [ -^rj exp ' " 

m=l 



ro I \ z \ ro 
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where = 2 [mW JiiUm )j ■ 

A/"' \ . . 

The function Jq ( 1 is strictly decreasing function if < r < ''o ■ Thus we can 

construct the partitions by the following partitions 

Bk,n = {{x,y)eC■^-^<^/^^T7<'-^^^^}, 0<i<n-l. 
{ n n } 

Now mes{Bk,n) = ~ g{-) where g{x) = nr'^x^, < a; < 1. 

Finally, the parameter of Poisson distribution is equal to 

1 

o = 2 ^rn^°^ Ji )^ 2nrQ J Joi^^i^ x)xdx = tt 



We used the following known relation aJo{a) = [aJi{a)]' [7,p.466] for calculation of 
the latter integral. 
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